Constraints on the axion-electron coupling for solar axions produced by Compton 

process and bremsstrahlung 
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The search for solar axions produced by Compton (7 + 6" — e~ + A) and bremsstrahlung-hke 
(e~ + ^ — >■ Z + + yl) processes has been performed. The axion flux in the both cases depends on 
the axion-electron coupling constant. The resonant excitation of low-lying nuclear level of ^^"Tm 
was looked for: A+^^^Tm -^^^'^Tm* ^^''^Tm +-f (8.41 keV). The Si(Li) detector and ^^''Tm target 
installed inside the low-background setup were used to detect 8.41 keV 7-rays. As a result, a new 
model independent restriction on the axion-electron and the axion-nucleon couplings was obtained: 
<?Ae X |(?Ajv + <?Aivl ^ 2.1 X 10^^*. In model of hadronic axion this restriction corresponds to the 
upper limit on the axion-electron coupling and on the axion mass qae x ttia < 3.1 x 10"^ eV (90% 
C.I.). The limits on axion mass are ttia < 105 eV and ttia < 1.3 keV for DFSZ- and KSVZ-axion 
models, correspondingly (90% c.L). 

PACS numbers: 14.80. Mz, 29.40.Mc, 26.65. -|-t 
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I. INTRODUCTION 

Axions arise as a result of the solution of the strong 
CP-problem proposed by Peccei and Quinn [1] . Wienberg 
and Wilczek showed that PQ-solution leads to the 
existence of a neutral spin-zero pseudoscalar particle with 
a non-zero mass. The axion mass as well as the strength 
of axion's coupling with ordinary matter is inversely pro- 
portional to the PQ-symmetry breaking scale Ja- The 
original PQWW-axion model with /a ~ (n/^Gf)"^/^ 
has exact predictions for axion coupling with photons, 
electrons and nucleons {qai, gAe and gAN )• The model 
has been excluded by the series of experiments with ra- 
dioactive sources, reactors and accelerators. 

Two types of the " invisible" axion models retained the 
axion in the form required for the solution of CP-violation 
problem, while suppressing its interaction with matter. 
These are the models of "hadronic" or Kim-Shifman- 
Vainstein-Zakharov (KSVZ) axion and the GUT or 
Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) axion 0,0]. 
The models of "invisible" axion have no certain predic- 
tions for /a, and as a result, no restrictions for the axion 
mass and the coupling constants. Moreover, the parame- 
ters of axion couplings are significantly model dependent. 

The hadronic axion does not interact with leptons and 
ordinary quarks at the tree level, which results in strong 
suppression of gAe constant through radiatively induced 
coupling l§\ . In some models axion-photon coupling may 
significantly differ from the original DFSZ or KSVZ cou- 
plings by a factor less then 10~^ Q. The effective axion- 
nucleon coupling depends on ratios of u-, d- and s-quark 
masses, axial pion-nucleon couplings F and D and poorly 
constrained flavor singlet coupling S. Moreover, the val- 
ues of gAN for the DFSZ axion depend on additional 
unknown parameter cos^ /? which is defined by the ratio 



of the Higgs vacuum expectation values @, Q . 

The axion mass niA (in eV units) in both models is 
given in terms of 7r° properties: 
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where /^r = 93 MeV is the pion decay constant; z = 
TTiu/mci = 0.56 and w = mu/ms = 0.029 are quark-mass 
ratios. The restrictions on the axion mass appear as a 
result of the restrictions on the coupling constants gA-y, 
gAe and gAN- 

If the axions do exist, then the Sun should be an in- 
tense source of these particles. Axions can be efficiently 
produced in the Sun by the inverse Primakoff conver- 
sion of the photons in the electromagnetic field of the 
plasma. The resulting axion flux depends on and 
can be detected by the Primakoff conversion of axions to 
photons in laboratory magnetic fields [l3| - [3 or by the 
coherent conversion to photons in the crystal detectors 
[13]- [2^. The expected count rate of photons depends 
on the axion-photon coupling as g\^. The upper limits 

on the gA-y are (10~^° — 10~®) GeV~^ for axions with 
mass less than 0.1 eV. 

There are two other possible mechanisms of axion pro- 
duction in the Sun: the reactions of solar cycle and the 
excitation of the low-lying energy levels of some nuclei by 
the high solar temperature. The attempts to observe the 
quasi-monochromatic axions emitted in nuclear magnetic 
transitions were performed in [13]- [l^. The reactions of 
the resonant excitation of nuclear levels, the axion-to- 
photon conversion and the axioelectric effect have been 
used for detection. 

In this letter we present the results of the search for 
axions emitted from Sun by the Compton process 7 + 
— )• + A and by bremsstrahlung + Z ^ + 
Z + A in the hot solar plasma. The cross sections of 
the both reactions depend on the axion-electron coupling 
constant g\^. The axions can be detected in the reaction 
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of the resonant absorption by -'^^^Tm nuclear target [ssj . 
The 8.41 keV 7-rays and conversion electrons produced 
by the de-excitation of the first nuclear level (Fig.l) can 
be registered. The detection probability of the axions is 
determined by the product x g\j^ which is preferable 
for small gAe or gA^y values. 

The results of laboratory searches for the axion as well 
as the astrophysical and cosmological axion bounds can 
be found in [io, |4l| . Constrains obtained with the solar 
axions remain of interest, even if they are less restric- 
tive than astrophysical arguments, because they are more 
comparable to the laboratory experiments. 



II. THE AXION SPECTRUM AND THE RATE 
OF SOLAR AXIONS ABSORPTION BY ^'^'^TM 
NUCLEUS 

If the axions or other axion-like pseudoscalar particles 
couple with electrons then they are emitted from Su n by 
the Compton process and by bremsstrahlung |42|-[47|. 
The expected spectrum of axions is calculated using the- 
oretica.1 predictions for the Compton cross section given 
in [H, [43 and the axion bremsstrahlung due to electron- 
nucleus collisions given in (s^l- The axion flux is deter- 
mined for radial distribution of the temperature T{r), 
density of electrons Ne{r) and nuclei Nz^a^t) given by 
BS05(OP) Standard Solar Model [13 based on high-Z 
abundances [53 |. 

The original photon flux for the Compton process is 
taken in accordance with Planck's law of black-body ra- 
diation; 



3.95 X 10^2 £:2 _ _ _ 
—^-j^^—LMY cm \ 



(2) 

The axion spectrum is found by integrating pho- 
ton spectrum dN^/dE^ and the Compton cross section 
da'^{Ey , Ea,tiia) / dEA over the radial distribution of T 
and Np-. 
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FIG. 1. 1,2 - the spectra of the axions produced by the 
Compton process and the bremsstrahlung, correspondingly 
{gAe ~ 10~^^, ruA =0). 3 - spectrum of the axions produced 
by Primakoff eflect (c/a^ = 10"^°GeV"^). The level scheme 
of ^®^Tm nucleus is shown in the inset. 



keV units. The corresponding solar axion luminosity is 
calculated to be La = g\e ^ 1-29 x lO^^L©, where L© 
is the solar photon luminosity. This value is in a good 
agreement with results obtained in [i^ [4^ . 

The spectrum of bremsstrahlung axions is calculated 
in the same way. The Maxwell-Boltzmann distribution 
for the kinetic energy of the scattering electrons at the 
temperature T{r) is used: 
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Since kT <^ nie we use the non-relativistic expression 
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The differential cross section for the axion 
bremsstrahlung process due to electron-nucleus col- 



for a{Ej) given in [48|, |49|. In this case there is a strong lisions da /dEA{Ee, Z) was calculated in [50|. The 



relation between axion and photon energy Ea — Ej and 
we omit the integration over Ej . The obtained spectrum 
at the Earth for gAe = 10^^^ and tua = is shown 
in Figdl line 1. In the (1 — 10) keV energy range the 
spectrum is parameterized with 1% accuracy by the ex- 
pression: 



d^A 
dEA 



cross section has complex form and it is proportional 
to Z^. In accordance with [34, 47] the contributions of 
electron-electron collisions to the axion bremsstrahlung 
emission are negligible at the some keV's axion energy. 
Therefore, only the number density of H, ''He, '^He, ^^C, 
i4]\j^ 16q g^j^jj 26pg nuclei in a given spherical shell of the 
solar interior at the radius r are taken into account: 
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where the value of the flux is given in cm '^s ^keV 
units and the value of the axion energy Ea is given in 
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where is cross section for the axion bremsstrahlung due 
to electron- nucleus collisions given in [HO] -The spectrum 
of solar bremsstrahlung axions calculated in assumption 
that gAe = 10~^^ and ttia = is given in FigUl hne 
2. The spectrum is softer then the Compton axion spec- 
trum, the maximum axion intensity corresponds to the 
0.6 keV energy and the average energy of axions is 1.6 
keV. The axion flux is well parameterized by the follow- 
ing expression (in cm~^s~^keV^^ units): 
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Due to the Compton process and bremsstrahlung 
the axion fluxes at 8.41 keV are d^A/dEA — g^e ^ 
1.13 X lO^^cm-^s-^keV"^ and d^A/dEA = g^e x 2.08 x 
lO'^^cni^'^s^^keV^^, correspondingly. For comparison, 
the spectrum of axions produced by the PrimakofF con- 
version of photons in electromagnetic field of the plasma 
is shown for gA^ = 10-i"GeV^^ (FigUl line 3). 

To take the dependence of axion spectra vs axion mass 
into account, the spectra were calculated for different val- 
ues of niA- The expected flux of 8.4 keV axions weakly 
depends on the possible values of ttia for uia < 5 keV, 
e.g. at TTiA — 5 keV the axion flux decreases by 30 %. 

As a pseudoscalar particle, the axion should be subject 
to resonant absorption and emission in the nuclear transi- 
tions of a magnetic type. In our experiment we have cho- 
sen the ^^^Tm nucleus as a target. The energy of the first 
nuclear level (3/2+) is equal to 8.41 keV, the total axion 
flux at this energy is 5^gXl.34xlO'^'^cm^^s~^keV~^. The 
8.41 keV nuclear level discharges through Afl-type tran- 
sition with i?2-transition admixture value of (5^=0.11% 
and the relative probability of 7-ray emission is 77 = 
3.79 X 10-3 [53. 

The cross-section for the resonant absorption of the 
axions with energy Ea is given by the expression that 
is similar to the one for 7-ray resonant absorption, but 
the ratio of the nuclear transition probability with the 
emission of an axion (uja) to the probability of magnetic 
type transition (w^) has to be taken into account. The 
rate of solar axion absorption by ^^^Tm nucleus will be 
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where (T07 is a maximum cross-section of 7-ray absorp- 
tion. The experimentally derived value of do-y for ^^^Tm 
nucleus is 0-07 = 2.56 x 10^^^ cm^ [H^]. A lifetime of 
the is^Tm first excited level is r = 5.89 ns fES], thus the 
width of energy level F = 1.13 x lO^^*^ keV. 

The ujA/i-^-y ratio calculated in the long- wave approxi- 
mation, has the following view [ssl . [56j : 
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Here, p^ and pa are the photon and axion momenta, 
respectively, ^0 = fip+^n ~ 0.88 and /Z3 = ^p— /z„ « 4.71 



are isoscalar and isovector nuclear magnetic momenta, /3 
and 77 are parameters depending on the particular nuclear 
matrix elements [IB, [ESl- In case of the ^^^Tm nucleus, 
which has odd number of nucleons and unpaired proton, 
in the one-particle approximation the values of /3 and rj 
can be estimated as /? « 1.0 and rj « 0.5. For the given 
parameters the branching ratio can be rewritten as: 
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In the KSVZ axion model the dimensionless isoscalar 
and isovector couplingconstants and g\j^ are related 
to fA by expressions [l, Q: 
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where M^r « 939 MeV is the nucleon mass. Axial- 
coupling parameters F and D are obtained from hyperon 
semi-leptonic decays with high precision: F=0.462 ± 
0.011, D= 0.808 ± 0.006 The parameter S charac- 
terizing the fiavor singlet coupling still remains a poorly 
constrained one. The boundaries (0.37 < < 0.53) and 
(0.15 < S < 0.5) were found in [5§| and ji^, accord- 
ingly. As a result the value of the sum (g^^v + 9\n) is 
determined within a factor of two, but the ratio w^/o;^ 
does not vanish for any value of parameter S. With 
S = 0.5 the numerical values of axion-nucleon couplings 
are: g\j^ = -4.03x lO-^m^ and gj^^ = -2.75x lO^^m^, 
mA is given in eV units. 

The values of g\j^ and ^^4^ for the DFSZ axion de- 
pend on the PQ charges of the u and d quarks [1, [§]. 
The charges X„ and have positive-definite values con- 
strained by relation Xu + Xd = 1. In case of ^^^Tm 
Ml-transition the value of (ujAj'-^'^y^^^'^ ratio lies within 
the interval - (0.11-=-1.74)(a;^/a;T,)^^^^ {S = 0.5). The 
lower and upper bounds of this interval are defined by 
values X^ = 0, AT^ = 1 and X„ = 1, = respectively. 

In accordance with ([7|-(fTU|. the rate of axion absorp- 
tion by ^^^Tm nucleus ([8]) dependent only on the coupling 
constants is (the model- independent view): 

Ra = 1.55 X 105.gl(5% + .g^^,)2(pA/P7)':S"'. (13) 

Using the relations between g\^^ g\j^ and axion mass 
given by KSVZ model (fTTj) . the absorption rate can be 
presented as a function of gAe and axion mass wia'- 

Ra = 5.79 X IQ-^"" g\^m\{pAlp^f .^-^ (14) 

In DFSZ-axion models the parameter gAe is associated 
with mass of the electron to, so that 
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where /3 is an arbitrary angle. If one sets cos^ (3=1 
the axion-electron couphng is related to axion mass as 
gAe='2'.8xl0^^^mA, where niA is expressed in eV units. 

The hadronic axion has no tree-level coupling to the 
electron, but there is an induced axion-electron coupling 
at the one-loop level 0: 



3a^Nm f E ^ Ja 2 4 + z + 'W ^ A 
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where N and E are the model dependent coefficients of 
the electromagnetic and color anomalies, A « 1 GeV 
is the cutoff at the QCD confinement scale. The nu- 
merical value of gAe for E/N — 8/3, which is char- 
acteristic for GUT models, and for = 3 is gAe = 

6.6 X 10-15 In - 14.6^ tua, where is ex- 

pressed in eV units. The interaction strength of the 
hadronic axion with the electron is suppressed by a factor 
^ c? . Because the coupling of DFSZ and KSVZ axions 
with electrons is much weaker than with nucleons the 
search for effect proportional to gAN x gAe is preferable. 
The dependence of axion absorption rate vs axion mass 
can be found through the relations ([T5|) and (fT6|) . 

The amount of observed 7-rays that follow the axion 
absorption depends on the number of target nuclei, mea- 
surement time and detector efficiency, while the proba- 
bility of 8.4 keV peak observation is determined by the 
background level of the experimental setup. 



III. EXPERIMENTAL SETUP 

To search for quanta with an energy of 8.41keV, the 
planar Si(Li) detector with a sensitive area diameter of 
66 mm and a thickness of 5 mm was used. The detector 
was mounted on 5 cm thick copper plate that protected 
the detector from the external radioactivity. The de- 
tector and the holder were placed in a vacuum cryostat 
and cooled to liquid nitrogen temperatures. A Tm2 03 
target of 2 g mass was uniformly deposited on a plex- 
iglas substrate 70 mm in diameter at a distance of 1.5 
mm from the detector surface. External passive shielding 
composed of copper, iron and lead layers was adjusted to 
the cryostat and eliminated external radioactivity back- 
ground by a factor of about 500. 

The experimental setup was located on the ground sur- 
face. Events produced by cosmic rays and fast neutrons 
were registered by an active shielding consisting of five 
plastic scintillators 50 x 50 x 12 cm in size. The rate 
of 50 /iS veto signals was 600 counts/s, that lead to w 
3% dead time. The Si(Li) detector was sectionalized into 
nine separate sections in order to lower the capacities 
of individual detectors and subsequently increase the re- 
sulting energy resolution. Every section was equipped 
with a charge-sensitive preamplifier with resistive feed- 
back, a shaping amplifier, and a 12-step analog-to-digital 
converter. 18 spectra consisting of 4096 channels from 



each detector (in anti- and in coincidence with the veto 
signals) are obtained. 

Though the detector amplifications were virtually the 
same, the energy calibrations were performed for each de- 
tector independently. Standard calibration sources ^^Co 
and ^"^^Am were used. Energy resolution in the integral 
spectrum for a 14.4 keV 7-ray line was a = 0.63 keV. 
The high energy resolution and accurate knowledge of 
the energy scale are crucial to our experiment because 
the energies of the characteristic X-rays of thulium are 
close to 8.41 keV. The most intense L-lines in the case of 
vacancy on K-shell have the next energies and intensities: 
7.18 keV (8.1%, L„i), 8.10 keV (5.2%, L/31) and 8.47 keV 
(1.6 %, L^2). 

The sensitive volume and the area of the Si (Li) de- 
tector were measured using the X-ray and 7-lines of a 
standard ^^^^Am source. The detection efficiency for an 
energy of 8.41 keV was estimated by numerical M-C sim- 
ulation, taking into account the self- absorption of 7-rays 
by the target. The simulation results were checked on 
a ^^^Am source placed behind the Tm203 target. The 
total 7-ray detection efficiency at an energy of 8.41 keV 
was £ = (6.16 ±0.30)%. 



IV. RESULTS 

Measurements were made over 31.8 days of live time by 
two hour series to monitor the stability of the Si(Li)- and 
the scintillation detectors. The integral energy spectrum 
of Si(Li)-detectors measured in the anticoincidence with 
signals of the veto scintillators is shown on the inset in 
Figl21 One can clearly identify peaks related to X-rays of 
thulium (K„i = 50.74 keV and K„2 = 49.77 keV). There 
were no statistically significant peaks in the spectrum of 
events correlated with the veto signals. 

Fig. 2 shows the detailed energy spectrum within the 
(6-;-20) keV interval, where the 8.41 keV is expected. The 
maximum likelihood method was used to define the in- 
tensity of the 'axion peak'. As the likelihood function 
N{E), we took the sum of the exponential function for a 
smooth background and a Gaussian: 



N{E) 



b X exp(c£^) 
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(17) 

The energy resolution a and peak position E^ were fixed, 
while the intensity Sa and background parameters a, 5, 
and c were free during the fitting. The result of the fit 
corresponding to the minimum of is shown in Fig. 2 by 
solid line. A standard method was used to set the upper 
limit on the 8.41 keV peak intensity: was determined 
for different fixed values of Sa while all other parameters 
were free. Obtained probability function P(x^(S')) was 
normalized to unity for 5* > 0. The upper limit estimated 
in this manner was Sum = 217 at a 90% confidence level. 
The expected number of registered 8.41 keV 7-quanta 
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FIG. 2. The Si(Li)-detector energy spectra measured in the 
anticoincidence with the veto signaL Solid hne shows the 
fitting resuh in the 6-20 keV range corresponding to the 
minimum x^- The spectrum in the (4 — 60) keV region is 
shown in inset. 

for the detection rate defined by (IT^ is 

Sa = eiiNie9TmTRA = 4.0 x lO^^i?^ < 5,„„, (18) 

where TViggTm = 6.23 x lO^^ - the number of i^STm 
nuclei, T = 2.75 x 10^ s - time of measurement, e = 
6.16 X 10~^ - detection efficiency and rj = 3.79 x 10^'^ - 
internal conversion ratio [53| . 

The upper limit on axions absorption rate by ^^^Tm 
nucleus Ra < 5.43 x 10~^'^s~^ set by our experiment 
limits the possible values of coupling constants gAe, 9 an 
and axion mass m^. According to (fT5)) and and 
taking into account the approximate equality of the axion 
and 7-quantum momenta (pA/p-y)^ — 1 for itia < 2 keV 
we obtain (at 90% c.l.): 

5AeX|(50jv+5i7v)l<2.1xlO-" (19) 

9Ae X < 3.1 X 10~^ eV (20) 

The restriction (jl9p is a model independent one on ax- 
ion (or any other pseudoscalar particle) couplings with 
electron and nucleons. The result (|20p presented as a re- 
striction on the range of possible values of gAe and itia 
(the relations (fTT|) and (fT^ between gAN and tha are 
used) allows one to compare our result with results of 



FIG. 3. The limits on gAe coupling constant obtained by 1- 
present work, 2- reactor experiments and Borexino !6J-[6y], 
3- beam dump experiments i67j, jjg,] , 4- ortopositronium de- 
cay Hi], 5 - CoGeNT m, 6- CDMS IH, 7 - solar axion 
luminosity [i^ 8 - red giant [4]]]). The areas of excluded val- 
ues are located above the corresponding curves. The expected 
values of qab and niA in DFSZ and KSVZ axion models are 
also shown. 



other experiments restricting gAe (FigO. The limits on 
gAe X iriA for DFSZ axion lie in the range (0.33 — 1.32) 
of the restriction (PO)) . 

The relation ^IU\\ excludes the region of relatively large 
values of gAe and niA possible in KSVZ and DFSZ mod- 
els. The strongest limit on gAe — 9 x 10"^^ corresponds 
to niA — 5 keV. The hypothesis of keV-scale bosons as 
possible dark matter candidates has been considered in 
j4S . i61j . More stringent limits on gAe was found under as- 
sumption that axion luminosity does not exceeds 0.1 Lq 
[igt . The recent constraints on keV-mass pseudoscalar 
dark matter by CoGeNT ^ and CDMS [63] are more 
restrictive than the solar luminosity limit. Reactor, beam 
dump and positronium decay experiments constrain the 
MeV-scale region of axion masses j63 - [69| . Very strong 
restriction on gAe can be obtained for axions with mass 
in (1.1 — 5.4) MeV range from the positron flux near the 
Earths atmosphere surface [3^. 

The obtained constraints (|19p and ([^0]) are valid in as- 
sumption that axions escape without restraint from the 
Sun. Axions leaving the center of the Sun pass through 
the matter layer with ss 6.8 x 10'^^ electrons/cm^. As 
the result the absorption of the axions due to inverse 
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Using the limit (|T9)) and relations p6|) and ([T5|) one can 
obtain the limit on the mass of KSVZ axion - ttia < 1-3 
keV and DFSZ axion niA < 105 eV (cos^ /3=1) (90%c.l.). 

A search for solar axions by resonant absorption or 
through axioelectric effect was made in [23|-[33|. The 
strongest limit for the hadronic axion mass (ruA < 151 
eV) was made for 14.4 keV axions emitted in the Ml 
transition of a ^^Fe nucleus [s^ . A significant advantage 
of our experiment is that for the Ml transition of ^^^Tm 
(as opposed to ^^Fe), uja/uj^ ratio depends weakly on the 
actual values of S and z. 

In our case, the uncertainty in S allows the upper limit 
to change from qac x niA < 2.7 x 10"'' eV {S — 0.7) to 
gAe X mA < 3.6 X lO"'' eV {S = 0.3)(Fig|ll). The value of 
u- and d-quark-mass ratio z = 0.56 is generally accepted 
for axion papers, but it could vary in the range (0.35—0.6) 
[40| . This uncertainty changes the obtained constraints 
insignificantly: from gAe x itia < 2.4 x 10"'' eV {z — 
0.35) to gAe X mA < 3.2 x 10"^ eV {z = 0.6). 

The sensitivity of the experiment depends on the total 
efficiency of registration which is77X£«2xlO"''in 
our case. This value can be increased significantly by 
introducing the Tm target inside the sensitive volume of 
detector. 



FIG. 4. The limit on qac x tua (in 10"^ eV units) versus the 
value of S parameter {p—l, j7=0.5, z=0.56, solid line). The 
doted, dash-doted and dashed lines correspond to the values 
of /3, -q and z changed by ±10%, correspondingly. 



Compton process A + e~ — ^ e" + 7 is lower than 10%, if 
gAe ^ 1x10"^. The other process associated with axion- 
electron coupling is the axioelectric effect A + e + Z ^ 
e + Z. The cross section of axioelectric effect was calcu- 
lated in [4^, 153 ■ The cross section has a Z^ dependence 
and solar abundance of Fe constrains the sensitivity of 
experiment by value gAe < 8 x 10"^. 



V. CONCLUSION 

We searched for the resonant excitation of the first 
nuclear level of ^^^Tm (8.41 keV) by axions formed in- 
side the Sun due to Compton and bremsstrahlung pro- 
cess provided by axion-electron coupling. A sectional- 
ized Si(Li) detector installed inside a low background 
setup was used to register 8.41 keV 7-quanta. As a 
result, we obtained a new model independent upper 
limit on the axion-electron and axion-nucleon couplings: 
gAe X \g% + fiwl < 2-1 X 10"" (90% c.L). In model of 
hadronic axion this restriction corresponds to the upper 
limit on the axion-electron coupling and on axion mass 
gAe X mA < 3.1 X 10"^ eV (90% c.L). 
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